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ABSTRACT 


From the data of a given aircraft, a maneuver model technique was 
developed which represented each maneuver type by a set of normalized 
parameter time histories, each independent of the given aircraft type. 

With proper denormalizing, the maneuver model can be applied in statis- 
tically predicting the maneuver-induced fatigue loads on any aircraft type. 
Thus the maneuver model could form a part of the structural design cri- 
teria. This technique was used to predict peak fuselage, wing, horizontal 
tail, and vertical shear load distributions from 450 hours of F-105D eight- 
channel oscillograph data. Favorable correlation of the predicted load 
peaks with load peaks calculated in the conventional fashion demonstrated 
the feasibility of the maneuver model technique. In an independent effort, 
a computer program was developed to determine the feasibility of auto- 
matically recognizing and classifying the maneuvers in eight-channel data 
digitized on magnetic tape. Using a magnetic tape bearing data from four- 
teen F-105D flights, an evaluation test showed that the program recognized 
90 percent of all maneuvers and correctly classified 90 percent of the ma- 
neuvers recognized. This test proved that such a program would reduce 
manual editing to less than one-tenth of the effort otherwise required. 
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STATISTICALLY PREDICTING MANEUVER LOADS 


FROM EIGHT-CHANNEL FLIGHT DATA 


By Larry E. Clay and 
Heber L. Short 
Technology Incorporated 


SUMMARY 


The objective of the effort reported was to demonstrate the feasi- 
bility of predicting maneuver loads fatigue spectra by using a set of 
normalized parameter data broken down by maneuver type and a set of 
normalizing factors broken down by maneuver type and by flight condi- 
tion. If feasible, this procedure could be utilized in design criteria 
where the normalized parameter data by maneuver type would be speci- 
fied as a ’'maneuver model. 11 To predict a maneuver fatigue load spectra 
for a given airplane requires estimating the type and number of maneuvers 
performed during the life of the airplane and the normalizing factor dis- 
tributions for each maneuver type. These estimated quantities would then 
be combined with the specified normalized data to predict the maneuver 
loads. 


The reported effort used a technique based on the maneuver model 
concept to predict the fuselage, wing, vertical tail, and horizontal tail 
loads from about 450 in-flight hours of eight- channel oscillograph data 
recorded during the peacetime operation of F-105D airplanes. First, 
for each of the four loads, all of the data were processed to calculate 
load time histories whose peaks were termed "observed loads. " Next 
the oscillograms were reviewed to establish twenty-three distinct and 
recognizable maneuver types. Then, from the maneuvers comprising each 
type, nine parameters were normalized in amplitude and time to form 
normalized data distributions by maneuver type. This process also 
yielded normalizing factors. Finally, in a simulation of the maneuver 
model technique to predict loads, the normalized data distributions were 
denormalized and the resultant data were used to predict the load peak 
distributions on the fuselage, wing, vertical tail, and horizontal tail. 

The favorable comparison of the predicted and observed load peak dis- 
tributions demonstrated the feasibility of the maneuver model technique 
to predict structural loads with acceptable accuracy. However, yet to 
be tested is the assumption that the normalized data is independent of 


the airplane type. The validity of this assumption would have to be proved 
before the maneuver model technique could readily use such data to pre- 
dict the loads of other airplane types. 

Since the practical use of the maneuver model technique on a large- 
scale data reduction basis depends on the extent to which present data 
editing and processing can be automated, an independent effort sought to 
develop a computer program capable of automatically recognizing and 
classifying maneuvers in digital eight-channel data. To evaluate the re- 
sultant program, magnetic tapes were prepared with digitized data simu- 
lating the eight-channel recordings made during fourteen flights of the 
F-105D airplanes. Results showed that the program recognized 90 per- 
cent of all maneuvers and classified correctly 90 percent of the maneuvers 
recognized. The effectiveness of the program indicated that computer 
processing could reduce the manual editing to less than one-tenth of the 
effort otherwise required. 


INTRODUCTION 


From the flight data of current aircraft, this study sought to evolve 
a maneuver model technique capable of predicting a maneuver flight loads 
spectrum. Such a spectrum would be used in designing new aircraft and 
in estimating the fatigue life of existing aircraft. 


Application of Flight Loads Data to Design Criteria 

The design of modern high-performance aircraft is an extremely 
complex task involving a delicate balance between the low weights re- 
quired for maximum performance and the high structural strength nec- 
essary to withstand the predicted aerodynamic and inertia loads. To 
verify predicted loads and to provide realistic data for the projection of 
the loads to be encountered by future aircraft is the purpose of flight 
loads programs. 

The design criteria for in-flight fatigue (Reference 1) consists pri- 
marily of a fatigue spectrum of cycles of percent design limit loads de- 
rived from three-channel (VGH) flight loads programs and a power spec- 
tral density of vertical gust velocities. These criteria, however, are 
limited to symmetrical vertical loads on the primary lifting surface (the 
wing). As discussed in Reference 2, such criteria for the other types of 
loads on other parts of the aircraft may be derived from the data acquired 
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in the so-called eight-channel flight loads data program. This study, 
therefore, used the eight-channel data recorded during an F-105D flight 
loads program. The parameters comprising this data are listed in Table 1. 


TABLE 1 

EIGHT-CHANNEL RECORDED PARAMETERS 
a x - longitudinal c. g. acceleration (positive forward) 
ay - lateral c. g. acceleration (positive right) 
a z - vertical c. g. acceleration (positive up) 
p - roll angular velocity (positive right wing down) 
q - pitch angular velocity (positive nose up) 
r - yaw angular velocity (positive nose right) 

- dynamic pressure 
P a - static pressure 


Three other required parameters — p (roll angular acceleration), q 
(pitch angular acceleration), and r (yaw angular acceleration) — are 
derived by differentiating the recorded analog traces of p, q, and r. 

And the airspeed, altitude, and Mach number are derived from the dy- 
namic and static pressures. Substituting these data in the rigid-body 
equations of motion yields a set of structural loads for any component 
of the aircraft. 

Because of the expense and complexity of recording and processing 
eight-channel oscillograph data, only a limited amount of this type of 
data now exists. However, digital recorders (predominately magnetic 
tape) will soon be in service to permit recording eight-channel data in 
volume and in a form compatible with high-speed digital computers. 

Since the equation for calculating a structural load is different for each 
point on the aircraft, the computing technique should be capable of cal- 
culating the structural loads from the recorded and derived parameters 
for any specific load spectrum without having to reprocess all the recorded 
data. 
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Maneuver Model Concept 


Unique to the maneuver model technique is the individual treatment 
of the flight loads data for each maneuver type. This technique is based 
on the belief that each maneuver of a given type will have the same se- 
quence of loads on a particular component of the aircraft structure, the 
only difference among the maneuvers being the amplitude and the duration 
of the loads. Therefore, through some suitable method of normalizing 
the amplitude and the time of the loads, all maneuvers of the same type 
would have the same time history of maneuver loads. 

With this basic concept of the maneuver model, two general pro- 
cedures would be followed to achieve the study objective: (1 ) a normal- 

ization procedure to develop the maneuver model, and (2) a denormal- 
ization procedure to apply the maneuver model. In the first procedure, 
the maneuver model would be formed after establishing the normalizing 
definitions. Then after selecting critical time slices and denormalizing 
the normalized data, the second procedure would apply the maneuver 
model to calculate the aircraft loads. 

Figure 1-a illustrates the steps for the normalization procedure 
leading to the development of the maneuver model. These steps are 
as follows: 

(1) Separate the recorded eight-channel data into maneuver types. 

(2) Normalize in time and amplitude each parameter trace in all 
maneuvers of each type. 

(3) For each maneuver type, combine all normalized parameter 
traces to form an average normalized parameter time history 
and to determine the distributions about the average trace at 
selected time slices. 

The normalized data resulting from step 3 comprises the maneuver 
model. Since the normalization of the parameter amplitudes is intended 
to make the data independent of the aircraft type, the maneuver model is 
directly applicable to structural design criteria. 

Figure 1-b illustrates the steps for the denormalization procedure 
to adapt and then apply the maneuver model to calculate, or predict, the 
fatigue load spectra for a specific structural member in a particular air- 
craft type. As noted in this figure, the application of the maneuver model 
requires denormalizing the normalized parameter data. To effect the 
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denormalizing, however, normalization factor distributions similar to 
those obtained in forming the maneuver model must be estimated for the 
particular airplane. The steps in Figure 1-b are as follows: 

(1) On the basis of the intended service life and mission mix of 
the airplane, estimate the total number of maneuvers of each 
type to be performed during the airplane life. 

(2) For each maneuver type, estimate the number of maneuvers 
in each flight condition (combination of variables such as air- 
speed, altitude, and weight). The flight condition will deter- 
mine the coefficients for the loads equations. 

(3) For each maneuver type and flight condition, estimate the nor- 
malizing factor distributions for each parameter. Available 
eight-channel flight data on similar types of airplanes will 
facilitate this estimation. 

(4) For each flight condition, calculate the coefficients for the 
loads equations. 

(5) For each maneuver type and flight condition, denormalize the 
average normalized parameter time histories and calculate 
an average load time history. Using a suitable peak criteria, 
locate the peaks, or critical times, on the average load time 
history. All maneuvers of each maneuver type are assumed 
to have load peaks only at these critical times. 

(6) For each maneuver type, each flight condition, and each load 
peak time, denormalize each parameter distribution and cal- 
culate the probability that the load was above each load level. 

(7) From the number of maneuvers in each maneuver type and 
flight condition combination and from the load probabilities 
found in step 6, calculate the distribution of the load peak 
values . 

As resulting from step 7, the composite of the load peak distri- 
butions for all maneuver types constitutes the maneuver-induced fatigue 
load spectra. 

In this feasibility study, the F-105D data was used to form the 
maneuver model. Then because the normalized data in this model was 
presumed to be independent of the aircraft type, the model could have 
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been used with the normalizing factors and loads equations for any air- 
craft type to predict the loads for such an aircraft. However, since the 
F-105D data was still the only available eight-channel data, the maneuver 
model was used with the F-105D normalizing factor distributions and loads 
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Limitations of the Maneuver Model 


The present maneuver model is limited to conventional fixed-wing 
aircraft types until its applicability to other aircraft types may be deter- 
mined by analyzing the recorded data from these other aircraft. Other 
possible limitations stem from two simplifying assumptions; (1) the 
independence of the normalized parameters at each selected normalized 
time and the normalizing factors, and (2) the insignificance of structural 
elasticity effects on the loads resulting from maneuvers. 


Limitations of Other Methods to Calculate Structural Loads 

References 2 through 5 present three other methods of calcu- 
lating maneuver loads on structural components from recorded flight 
data. Each of these methods requires statistical techniques because of 
the huge data sample and the variation in the sample due to the uncon- 
trollable effects of different pilot techniques, atmospheric turbulence, 
geographic topology, and weather conditions. 

The method given in References 2 and 3 utilizes discrete samples 
of the parameters (usually the parameter peak values). This method has 
some serious limitations, two of which are described below; 

The first shortcoming of this method lies in its inability to yield 
the number of load cycles. Two variables essential to a fatigue load 
spectrum are the number of load cycles and the magnitude of the peaks 
in these cycles. Although this method can calculate a probability of 
exceeding given structural loads in a given number of flights from the 
recorded and derived parameters, it cannot determine the number of 
load cycles from them. Knowing the number of recorded parameter 
peaks does not suffice since the number of peak loads on some of the 
aircraft structural components does not necessarily correlate with it. 

For example, since vertical tail loads can be considered basically a 
function of the lateral load factor ny and the yaw angular acceleration r, 
it would seem reasonable to correlate the number of vertical tail load 
peaks with the number of peaks of each of these parameters. However, 
Figure 2, a time history for part of a left turn, shows an instance where 
ny and r peak at slightly different times with only a single corresponding 
vertical tail load peak. In other instances, a single vertical tail load 
peak still appears when only one of the ny and r parameters peaks. Con- 
sequently, the total number of vertical tail load peaks is probably more 
than the total number of either the ny or r peaks but less than the sum 
of the totals . 
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Figure 2. —Sample time histories of tail load, yaw acceleration, 
and lateral acceleration 


The second shortcoming of this method derives from the inherent 
assumption, not always true, that apeak load occurs when one of the 
basic parameters peaks since the method samples all parameters only 
at the instant that one of them peaks. However, as shown in Figure 2, 
the calculated vertical tail load peak occurs between the n^ and r peaks. 
During the processing of the data, the ny and r values corresponding to 
the vertical tail load peaks are discarded unless the load peak occurs 
simultaneously with a parameter peak. Thus the magnitude of loads 
derived from the processed data is not necessarily correct. 

References 4 and 5 present, respectively, the Monte Carlo and 
power spectral analysis techniques to calculate the maneuver loads. The 
shortcoming of these techniques is due to their random selection of meas- 
ured or calculated parameter values to calculate the loads. Since. a pilot 
performs a maneuver to change the aircraft attitude or to meet his vari- 
ous mission requirements, like maneuvers will often be repeated in 
succession. With such controlled operation the resultant parameter 
peaks are certainly not random with respect to time. 

The shortcomings of these other methods led to the concept of the 
maneuver model technique with its unique individual treatment of ma- 
neuver types to provide a more reasonable method of calculating maneu- 
ver loads on aircraft structural components. 
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Preliminary Feasibility Study of the Maneuver Model Technique 

As reported in Reference 6, a preliminary feasibility study of the 
maneuver model technique was conducted under Contract NASw-970. 
After the F-105D eight-channel data were processed for a single maneu- 
ver type, the technique was used to calculate the loads on the wing, the 
horizontal tail, and the vertical tail. Then these predicted loads were 
compared with the observed loads derived from the time histories of the 
recorded parameters. The favorable results led to the larger scale 
feasibility program reported here. 
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SYMBOLS 

aircraft 

correction constant for the ith normalized time slice 

2 

longitudinal acceleration, feet/second — positive forward 

2 

lateral acceleration, feet/second — positive right 
normal acceleration, feet/ second — positive up 
wing span, feet 

mean aerodynamic chord, feet 
center of gravity 

coefficients of the loads equations 

lift coefficient per degree of a 

lift coefficient per degree of i^ 

lift coefficient on the wing and fuselage at a = 0 

rolling moment coefficient per degree of 3 

rolling moment coefficient per degree of 6 ^ 

rolling moment coefficient per degree of 

rolling moment coefficient due to 6g per degree of 6 a 


'*v 


c h 


J m r 


m i*p 

mq 


m r 




°HT 


'nc 


n &A 


j n T 


' nr 


' y P 

"ypvT 

: Y6 A 


'Y6 


R 


Y r 

D/S 


D/S/S 


rolling moment coefficient per radian of pb/2V 

rolling moment coefficient per radian of rb/2V 

pitching moment coefficient per degree of a 

pitching moment coefficient per degree of i'p 

pitching moment coefficient per radian of qc72V 

pitching moment coefficient at a = 0 

horizontal tail pitching moment coefficient at a = 0 

yawing moment coefficient per degree of 3 

yawing moment coefficient per degree of 6 A 

yawing moment coefficient per degree of 

yawing moment coefficient per radian of pb/2V 

yawing moment coefficient per radian of rb/2V 

side force coefficient per degree of 3 

vertical tail side force coefficient per degree of 3 

side force coefficient per degree of 6 A 

side force coefficient per degree of 

side force coefficient per radian of rb/2V 

degrees/ second 

degrees/ second^ 

fraction of a at horizontal tail a 

2 

gravitational constant, feet/second 
horizontal tail incidence, degrees 

2 

moment of inertia about the x-axis, slug-feet 


moment of inertia about the y-axis, slug-feet 


angular roll rate, degrees / second— positive right wing down 

2 


I_ moment of inertia about the z-axis, slug-feet 11 

z 

2 

I xz product of inertia, slug-feet 

n x , N x ,NX: longitudinal load factor (a x /g) 

n y , N y ,NY lateral load factor (a y /g) 

n z , N z , NZ normal load factor (a z / g + 1.0) 

An , DNZ incremental normal load factor (n z - 1.0) 

Z 

p.p 

p, PDOT angular roll acceleration, degrees/second 
static pressure, inches of mercury 
dynamic pressure, inches of mercury 
probability of indicated event 

angular pitch rate, degrees/second — positive nose up 

q, QDOT angular pitch acceleration, degrees/second 

r, R angular yaw rate, degrees/ second — positive nose right 

2 

r, RDOT angular yaw acceleration, degrees/second 
wing area, feet^ 
time, seconds 
ith normalized time slice 
true airspeed, feet/second 
equivalent airspeed, knots 
fuselage shear load, pounds — positive up 
maneuver fuselage shear load (Vp - Vp^ _ P oun ds 


r a 

Pd 

Pr[ ] 

q, Q 


s 

t 

^i 

V 

V e 

V F 

AVx 
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V 


HT 


AV 


HT 


v VT 
AV 


VT 


W 


w 


AV 

w 

w 6 

x 

x 

x i 

( x i>< 


x j 


a 

3 

6a 

6 r 

6 S 


horizontal tail shear load, pounds — positive up 

maneuver horizontal tail shear load (Vpj^p - _ ^), pounds 

vertical tail shear lead, pounds — positive right 

maneuver vertical tail shear load (Vyp - Vy»j^ - P oun< ^ s 

wing shear load, pounds — positive up 

maneuver wing shear load (Vw - V w ), pounds 

(n z = 1) 

aircraft instantaneous weight, pounds 
weight on right outboard wing pylon, pounds 

longitudinal coordinate — positive forward 

average of the j peak normalized parameter values 

normalized parameter value at ith normalized time slice 

corrected normalized parameter value at ith normalized 
time slice 

average normalized parameter value at ith normalized time slice 

peak normalized parameter value for the jth maneuver 

lateral coordinate — positive right 

normal coordinate — positive down 

angle of attack, degrees — positive nose up 

angle of sideslip, degrees — positive nose right 

aileron deflection, degrees — positive right aileron up 

rudder deflection, degrees — positive left 

spoiler deflection, degrees — positive right spoiler up 

local atmospheric density, slugs /foot'* 
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3 

n standard sea level atmospheric density, slugs /foot 

K o 

a density ratio (p/Po) 

Subscripts 

c this symbol indicates corrected values in normalized data 

i, j dummy indices 

(n z = 1) this symbol indicates level flight values for loads 

P this symbol indicates approximate expressions for loads with 

second-order terms eliminated 


SECTION I 
PROCEDURES 


A. Type and Source of Data Used 

Eight-channel data collected with oscillograph recorders in F-105D 
aircraft were chosen to demonstrate the application of the maneuver model 
techniques. These data represent 450 hours of the normal peacetime 
flight of F-105D's operated from the following bases: Bitburg Air Base, 

Germany, Wheelus Air Base, Libya; Kadena Air Base, Okinawa; and 
Nellis Air Force Base, Nevada. Of this data, 250 hours were previously 
reduced and reported in Reference 7. 

The recording system included the following sensors: three strain 

gage accelerometers, each mounted at the airplane's center of gravity 
and aligned with one of the three major axes of the airplane, to measure 
the normal (vertical), lateral, and longitudinal acceleration; three poten- 
tiometer rate gyros, each aligned with one of the three major axes of the 
aircraft, to measure the angular rates around these axes, that is, pitch, 
roll, and yaw; and two pressure transducers, connected to the airplane's 
pitot-static system, to measure the static and the dynamic pressure. 
According to the manufacturers of these instruments, the accuracies for 
the accelerometers and the pressure transducers were ± 1 percent of 
full scale, and those for the angular rate gyros were ± 3 percent of full 
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scale. Only flight data were recorded since the oscillograph was not 
turned on until after takeoff when the landing gear was retracted and 
turned off before landing when the landing gear was extended. 


B. Calculations and Definitions 

1. Parameters and variables. —To facilitate the descriptions 
given in this report, the terms "parameters" and "flight condition 
variables" are defined as follows: "parameters" denotes the airplane 

motion variables n x , n-y, n z , p, q, r, p, q, and r. And "flight condition 
variables" designates any or all of the variables defining a given set of 
loads equations; some of the flight condition variables are weight, mo- 
ments of inertia, airspeed, altitude, Mach number, stores configuration, 
and fuel weight. The term "flight condition" is used to identify any group 
ing of flight condition variable values into the complete set required to 
define the loads equation coefficients. 


All flights were categorized by five mission types, as listed in 
Table 2; and each flight was divided into mission segments, according 
to the fourteen classifications given in Table 3. The flights were also 
categorized according to the arrangement of the fuel tanks and the stores 
on the outboard wing pylons. For this purpose, one configuration for the 
external fuel tanks and centerline stores and a second for the stores on 
the outboard pylon were arbitrarily chosen. Table 4 defines these store 
configurations and gives the identifying codes for them. 


TABLE 2 

F-105D MISSION TYPES 

(1) Special Weapons Delivery 

(2) Conventional Bombing and Ground Gunnery 

(3) Air-to-Air Gunnery 

(4) Air Tactics and Test Flights 

(5) Instruments and Navigation 
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TABLE 3 


F-105D MISSION SEGMENT TYPES 


(1) 

Ascent 

(8) 

Rockets 

(2) 

Cruise-out 

(9) 

Air-to-Air Gunnery 

(3) 

Loiter 

(10) 

Air-to-Ground Gunnery 

(4) 

Low-Angle Bombing 

(ID 

Refueling 

(5) 

High-Angle Bombing 

(12) 

Training 

(6) 

GAM Delivery 

(13) 

Clean Cruise 

(7) 

Special Weapons Delivery 

(14) 

Descent 

Aire 

raft gross weight at liftoff. 

when the recorder was turned 


and just before landing, when the recorder was turned off, was based on 
the flight-line-departure and flight-line-return fuels logged by the field 
technicians. For the liftoff weight, 2300 pounds was subtracted from 
the flight-line-departure weight to adjust for the fuel burned during start- 
up, taxi, and takeoff. And for the "touchdown" weight, 400 pounds was 
added to the flight-line-return weight to account for the fuel consumed 
during landing, taxiing, and parking. With the fuel consumption rate 
during flight assumed constant, a rate of fuel-weight loss was then com- 
puted to find the instantaneous gross weight of the aircraft. When ex- 
ternal stores were dropped or the aircraft was refueled in flight, the 
weight was adjusted accordingly. 

As supplied by Reference 8, the moments of inertia for a single 
configuration and a specific gross weight served as base values. As the 
stores configuration and the fuel weight varied from the base configu- 
ration and weight, estimated increments were added to or subtracted 
from the base moments of inertia. The sequence of fuel tank usage in 
each flight was assumed to be that recommended in Reference 9- When- 
ever external stores were dropped in flight, the moments of inertia were 
modified to the new configuration. 

After the oscillogram trace displacements were digitized, they 
were converted by linear calibration slopes to the corresponding physical 
units, that is, units of g's, degrees per second, inches of mercury, and 
minutes. Calibration pulses recorded prior to each flight provided the 
means of adjusting the appropriate calibration slopes. 

The 19 59 Standard Atmosphere tables were used to derive the indi- 
cated pressure altitude, calibrated airspeed, and indicated Mach number 
from the dynamic and the static pressure. Then the conversion tables in 
Reference 10 along with the corrections for the pitot- static position error 
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in Reference 9 were used to convert these parameters to pressure alti- 
tude, equivalent airspeed, and true Mach number. 


From each of the three angular rate traces of p, q. and r, the cor- 
responding angular accelerations of p, q, or r were derived by taking the 
time derivative of the mid-point of the parabola formed through each set 
of three consecutive readings, where the first reading in each set was the 
second in the preceding set. In other words, at time t = 2, 


P2 


?3 ~ Pi 
*3 " fc l 


TABLE 4 

STORE CONFIGURATION CODES 
External Fuel Tank and <k Store Configuration 

Code 

1 clean — no auxiliary tanks or <£, stores 

2 BDU on <£, only 

3 bomb -bay tank only 

4 bomb -bay tank and 450 -gal. tank on each inboard pylon 

5 450- gal. tank on each inboard pylon 

6 BDU on <&, and 450-gal. tank on each inboard pylon 

7 bomb-bay tank, 650-gal. tank on <& pylon and 450-gal. tank on each 
inboard pylon 

8 BDU, 650-gal. tank on (£, pylon and 450-gal. tank on each inboard pylon 

9 multiple stores ejector in bomb bay and 450-gal. tank on each inboard 
pylon 

N. B. : BDU refers to any store of about 2000 lb. 


Outboard Pylon Store Configuration 


Code Left Outboard Pylon 

1 none to 150 lb. 

2 150 to 500 lb. 

3 none to 150 lb. 

4 150 to 500 lb. 

5 above 500 lb. 

6 none to 150 lb . 

7 none to 150 lb. 

8 above 500 lb. 


Right Outboard Pylon 
none 

50 to 150 lb. 

50 to 150 lb. 

150 to 500 lb. 
none 

150 to 500 lb. 
above 500 lb. 
above 500 lb. 


17 


2. Loads equations. —To calculate load values from the measured 
eight-channel data at 1/5-second intervals, a set of loads equations were 
developed and thresholds for the load peaks were defined. Hereafter the 
calculated load values are referred to as "observed loads. " Appendix F 
explains the development of the loads equations and gives the observed 
load peak definitions. 

3, Maneuver types . — To distinguish types in recorded data re- 
quires identifying the characteristics of each as discussed in the time 
histories of the parameters. In the study of the F-105D data to evolve 
distinguishable maneuver types, the criterion was the manner in which 
the airspeed and altitude traces increased or decreased and the other 
traces deflected from their normal positions during straight and level 
flight. As a result, twenty-three maneuver types were judged distinct 
and identifiable as listed in Table 5. 


TABLE 5 

F-105D MANEUVER TYPES 


(1) Descending left turn 

(2) Level left turn 

(3) Ascending left turn 

(4) Descending right turn 

(5) Level right turn 

(6) Ascending right turn 

(7) Symmetrical pull-up 

(8) Right rolling pull-up 

(9) Left rolling pull-up 

(10) Right roll 

(11) Left roll 

(12) Longitudinal deceleration 


(13) Longitudinal acceleration 

(14) Left yaw 

(15) Right yaw 

(16) Left wing rock 

(17) Right wing rock 

(18) Left cloverleaf 

(19) Right cloverleaf 

(20) Symmetrical pitch-down 

(2 1 ) Inside loop 

(22) --- (23) 

(24) Left four-point roll 
(2 5) Right four -point roll 
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The absence of maneuver types for Nos. 22 and 23 is explained as 
follows: Since the evolvement of the maneuver types was concurrent with 

the development of the computer program, an optimistic twenty-five num- 
bers were allocated to ensure the inclusion of all maneuver types dis- 
cerned. When the two four-point roll maneuvers were established, they 
were judged sufficiently different from the others then determined to be 
numbered last. Then after all maneuver types were established, they 
were not renumbered to Nos. 22 and 23 because of the effort required 
to change the computer program and the data already so identified. 
Appendix B describes and illustrates the maneuver types at length. 


4. Normalization procedure to develop maneuver model. — 

a. Normalization definitions: A set of criteria was arbitrarily 

established to normalize the data for each maneuver type. This normal- 
izing was intended so that the time histories of the parameters for each 
maneuver type could be described independently of the airplane type, 
the magnitude of the parameter deflections, and the abruptness of the 
maneuver. To normalize the time of each maneuver, two easily recog- 
nized anchor points were chosen as outlined in Table 6 and were assigned 
two normalized time values. Then, on the basis of these two anchor 
points, the original time scale was linearly transformed into the normal- 
ized time scale. Next the amplitude of each parameter trace in each ma- 
neuver was normalized independently by dividing each reading by the 
maximum deflection of the trace. (Note that in normalizing the n z trace, 
the normalized data and maximum deflection were calculated from the 
incremental normal load factor An z . ) These maximum deflections are the 
"normalizing factors" referred to in the introduction and henceforth will 
be called "maximum absolute parameter values. " These maximum abso- 
lute parameter values were stored for each maneuver. At the end of 
the data processing, the maximum absolute parameter values formed a 
set of nine maximum absolute parameter distributions for each maneu- 
ver type. 


TABLE 6 

TIME-NORMALIZING CRITERIA 



First Anchor Point 


Second Anchor Point 


Man. 

Type 

Description 

Normalized 

Time 

Description 

Normalized 
T ime 

1-3 

min. p prior to max. n z 

0 . 1 

max. p in last . 4 of maneuver 

0. 9 

4-6 

max. p prior to max. n z 

0 . 1 

min. p in last . 4 of maneuver 

0. 9 
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TABLE 6. —Continued 


TIME-NORMALIZING CRITERIA 



First Anchor Point 


Second Anchor Point 


Man. 

Type 

7 

Description 

Normalized 

Time 

Description 

Normalized 

Time 

max. q at or prior to max. n z 

0. 2 

(max. An z )/2 after max n z 

0.45 * 

8 

max. q at or prior to max. n z 

0. 1 

min. p in last . 3 of maneuver 

0. 9 

9 

max. q at or prior to max. n z 

0. 1 

max. p in last 3 of maneuver 

0.9 

10 

(max. pl/2 before . 6 of maneuver 
and before max. p 

0.2 

(max. p)/2 after 5 of maneuver and 
after max. p 

0. 75 

11 

(min. p}/2 before . 6 of maneuver 
and before min. p 

0.2 

(min. pl/2 after . 5 of maneuver and 
after min. p 

0. 75 

12 

first n x decrease of . 015 or greater 
in first . 4 of maneuver 

0. 1 

first increase of . 035 or greater 

after min. n x in last . 5 of maneuver 

0. 75 

13 

max. n x in first . 1 5 of maneuver 
within 3 readings of first n x 
increase of 100% or greater 

0. 05 

min. n x in last . 4 of maneuver 
within 3 readings of first n x 
decrease of 75% or greater 

0. 95 

14 

first relative min. r in 
first . 5 of maneuver 

0. 15 

first relative max. r following first 
anchor point in first 75 of maneuver 

0 4 

15 

first relative max. r in 
first . 5 of maneuver 

0. 15 

first relative min. r following first 
anchor point in first . 75 of maneuver 

0 4 

16 

first relative min. p with 
a value below -15°/sec in first 
35 readings of maneuver 

0. 15 

first relative max. p following 
first anchor point with a value 
above + 1 5°/ sec . 

0. 4 

17 

first relative max. p with 
a value above + 15°/sec in first 
35 readings of maneuver 

0. 15 

first relative min. p following 
first anchor point with a value 
below -15°/ sec 

0. 4 

18 

min. p in first . 1 5 of maneuver 

0. 05 

max. n z in last . 5 of maneuver 

0. 85 

19 

max. p in first , 15 of maneuver 

0. 05 

max. n z in last . 5 of maneuver 

0. 85 

20 

min. q in first . 5 of maneuver 

0. 3 

(min. An z )/2 after min. n z 

0. 7 

21 

max. n z in first , 5 of maneuver 

0. 3 

max. n z in last . 5 of maneuver 

0. 8 

24 

first relative min. p with a value 
below -25°/ sec in the first . 5 
of maneuver 

0. 15 

fourth relative min. p with a value 
below -25° /sec and separated from 
first anchor point by three relative 
max. p's — each at least 20% greater 
than its preceding minimum 

0. 65 

25 

first relative max. p with a value 
above +25°/ sec in the first . 5 
of maneuver 

0. 15 

fourth relative max. p with a value 
above +25° /sec and separated from 
first anchor point by three relative 
min. p's — each at least 20% less 
than its preceding maximum 

0. 65 
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Each maneuver type has a definite, uniform pattern of parameter 
deflections which reflect the change in aircraft controls to effect the 
desired maneuver. In the midst of this pattern,' however, are some ran- 
dom low-magnitude deflections caused by atmospheric turbulence and 
inadvertent control changes. Whenever all parameter deflections in a 
maneuver are so low that the random low-magnitude deflections appear 
significant, the normalizing of that maneuver would yield meaningless 
normalized data. The following thresholds, therefore, were established 
to determine whether the parameters defining the characteristic pattern 
of a maneuver type were sufficiently large to warrant normalizing the 
maneuver. 


Maneuver Type 

Turns (Types 1,2,3, 4, 5, 6, 18, and 19) 

Pull-ups (Types 7,8,9, and 21) 

Rolls (Types 10,11,16,17,24, and 25) 
Accel, and Decel. (Types 12 and 13) 
Yaws (Types 14 and 15) 

Pitch-down (Type 20) 


Thresholds 


n z - 2- 0 g 

-5 g r s 5°/sec 
- 3. 0 g 

n z - 2 ■ 0 g 

- • 2 - n x § . 2 g 

-5 ^ r 5 5°/sec 
no threshold 


- . 1 ^ n y g . 1 g 
-30 g p s 30°/ sec 


-30 S p § 30° /sec 


- . 1 ? ny s . i g 


After the maneuvers of Types 1 through 7 and 10 through 19 had 
been normalized in some 200 hours of data to yield an adequate statis- 
tical sample of these types, no more like maneuvers were normalized. 

The normalized parameter distributions should be formed at nor- 
malized times corresponding to the maximum and minimum loads. But 
since these loads depend on the aircraft configuration, the structural 
point of interest, and the relative magnitude of several parameters, the 
normalized times of their occurrence cannot be determined in advance. 
Therefore, to provide normalized parameter distributions at sufficient 
times to permit correlating them later with the time of any load, the 
normalized parameter distributions at several normalized time slices 
were formed. A preliminary review of the number and duration of the 
load peaks in the F-105D data indicated that twenty-five normalized 
time slices located near expected load peak times would suffice. 
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b. Formation of maneuver model: The statistical model to pre- 

dict maneuver loads categorizes maneuvers according to 23 types. As 
explained in the previous section, the maneuvers of each type found in 
the sample of F-105D eight-channel data were normalized both in time 
and in amplitude. Then for each of 25 normalized time slices, the nor- 
malized data for each of the nine parameters were stored in frequency 
tables. Similarly, the maximum absolute values of each parameter for 
each maneuver were stored in frequency tables along with the associated 
maneuver type and flight condition. 

The selection of a flight condition and maneuver type determined 
the exact maximum absolute parameter frequency distributions needed 
for loads predictions. These frequency distributions were converted into 
probability distributions by dividing each frequency by the number of ma- 
neuvers in the data sample of the given flight condition and maneuver 
type. Then later this same number of maneuvers was used in the pro- 
cess of converting the probabilities of loads into predicted load fre- 
quencies. Similarly, the frequencies stored in the normalized parame- 
ter data tables were converted into probabilities by dividing by the 
number of maneuvers which were normalized for the maneuver type. 

The normalized data consists of (1) a set of average normalized 
time histories of the nine aircraft parameters (n x , n y , n z , p, q, r, p, 
q, and r) for each maneuver type, and (2) a description of how the nor- 
malized parameter values would be expected to vary from the average 
at each of the 25 normalized time slices. An average normalized 
parameter trace was the time history formed by 25 values, each equal 
to the mean of the distributions of the normalized parameter values at 
one of the 25 normalized time slices. For each maneuver type, nine 
average normalized parameter time histories (one for each parameter) 
were formed. Maximum positive and maximum negative values of an 
average normalized parameter time history are the maximum and mini- 
mum of the 25 means, respectively. 

Each time an individual maneuver was normalized, it contributed 
a single frequency count to each of the 25 normalized distributions at the 
25 normalized time slices for each of the nine parameters. The nor- 
malizing procedure was designed to align all peaks of each parameter 
and the time slices were chosen so that some of the normalized parame- 
ter distributions would be oriented as close as possible to all peaks. 
However, because of the variation inherent in the data, not all parame- 
ter traces peaked precisely at a chosen time slice. Consequently, the 
maximum average normalized parameter value was the mean of a dis- 
tribution where not all values were peaks. Obviously, such a maximum 
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average normalized parameter value would have heen higher if all pa- 
rameter traces had peaked at that time slice. 

The nine average normalized parameter time histories for each 
maneuver type were used to calculate an average load time history for 
that maneuver type. Then, all load peaks were assumed to occur at the 
same normalized times as those in the average load time history. In 
reality, however, it would be expected that the individual load peaks 
would be misaligned in the normalized time scale as much as the nor- 
malized parameter peaks were misaligned. Thus the load distributions 
calculated at the normalized time slice corresponding to the average load 
time history peaks would be an underestimation of the true load peak dis- 
tribution by an amount proportional to the underestimation of the pertinent 
maximum average normalized parameter values. To compensate for 
this, each parameter distribution at each of the 25 normalized time slices 
was adjusted. 

The procedure for adjusting the normalized parameter values is 
illustrated by the exaggerated example in Figure 3. This figure shows 
the superimposed normalized traces of the same parameter recorded in 
three maneuvers and four of the normalized time slices. Each trace has 
a normalized peak value £ j . Each normalized time slice ti has an average 
normalized parameter value which is the mean of the distribution of the 
three x^ values at that time slice. The maximum average normalized 
parameter value is X£. The adjustment procedure was designed to correct 
this maximum average value to the average x of the j individual normal- 
ized trace peak values defined by 


— r-> A ^ A A 

A Zx i xi + X2 + x 3 

n 


where 


j = 1, 2, ... n 

n = 3 = number of normalized maneuvers 

A positive correction constant A£ at normalized time slice t2 (where 
the maximum average normalized parameter value occurred) was calcu- 
lated as follows: 


A - ,£&j - ' 

A 2 = x - x 2 = - x 2 ) 
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NORMALIZED x 



NORMALIZED t 

Figure 3. —Three exaggerated normalized traces of the same 

parameter to illustrate correction for underestimated 
load peak distributions 
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And correction constants A[ for the other normalized time slices t^ (with 
positive average normalized parameter values) were defined by 


Ai = 


Xj_ 


x 2 




Then every element x^ of the normalized parameter distribution at the 
normalized time slices tj was corrected by adding the correction con- 
stant Ai as follows: 

( x i)c “ x i + ^i 


where 


( x i) c = element of corrected distribution. 

Negative correction constants defined for each normalized time 
slice t^ where the average normalized parameter value 3q was negative 
were calculated in the same manner except the xj peaks were the mini- 
mums of the individual normalized parameter traces and the t^ normal- 
ized time slice corresponded to the minimum average normalized pa- 
rameter value x 2» 


5. Denormalization procedure to apply maneuver model. — 

a. Selection of critical time slices (load peaks ): For each pa- 

rameter in a given maneuver type, an average denormalized time history 
of the parameter was calculated by finding at each of the 25 time slices 
the average corrected normalized parameter value from the corrected 
normalized distribution. Then the average value from the maximum ab- 
solute parameter distribution was formed. The product of these two 
averages gave the average denormalized value of a parameter at each 
time slice. These average parameter values at each of the 25 time slices 
were used in the loads equations given in Appendix A to compute the aver- 
age load time history for each load type. 

The coefficients in the loads equations are functions of the flight con- 
dition variables. The value assigned to a given flight condition variable 
was found by calculating a weighted average of the mid-values in 20 ranges 
of that flight condition variable. The weight used for each mid-value was 
the number of maneuvers which fell in the corresponding range of the 
variable. 
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The calculated average load time histories for each maneuver type 
were delta loads. The steady load values were constant for each flight 
condition because they are independent of any parameter variations. The 
average delta load time histories were necessary to define the predicted 
load peak criteria since the observed load peak criteria were based on 
delta load time histories. Throughout the load prediction calculations, 
the delta load peak distributions were predicted and the corresponding 
ranges of steady loads were carried along. 

Since it was assumed that the average load history describes the load 
history of all maneuvers of this type, all of the load distributions were 
predicted at normalized time slices which were peaks of this average load 
history. The choice of which normalized time slices and of the threshold 
values to be associated with the time slices was determined by the peak 
criteria placed on the observed peak loads. 

Remembering the threshold values placed on observed peak loads, 
criteria for choosing the critical normalized times of predicted peak loads 
from the average load trace were established as follows. First, every 
occurrence of a maneuver type was assumed to have a load trace with the 
shape (not the magnitude* of the average load trace. Thus, the first part 
of the observed peak criteria — the trace must have a rise and fall of 50 
percent of the peak value or greater — could be applied to the average load 
because the shape of the trace was independent of the magnitude. 

Thus every point on the average load trace preceded by a rise of at least 
50 percent of its value and followed by a fall of at least 50 percent of its 
value became a tentative positive load peak. (Tentative negative load 
peaks had the same criteria except for the sequence of a fall and then a 
rise. ) The rest of the observed peak criteria — the peak value must be 
outside threshold and the amount of rise and fall must at least equal the 
threshold value — required that the magnitude of the load trace be known. 
The tentative peaks were divided into two groups: (1) tentative primary 

load peaks where the rise and fall of the trace were each equal to or 
greater than the peak value, and (2) tentative secondary load peaks where 
the rise or fall or both were less than the peak value but both were greater 
than 50 percent of the peak value. For the tentative primary load peaks, 
a threshold value was established which was equal to the observed load 
peak thresholds. Later during the load probability calculations, when- 
ever a tentative primary load peak time slice had a predicted load fre- 
quency distribution with load values equal to or greater than this threshold 
value, it satisfied the peak criteria. For the tentative secondary peaks, 
a different threshold value was calculated so that when the peak value was 
equal to it, both the rise and fall of that peak were at least equal to the 
observed load peak threshold. The secondary load peak thresholds were 
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at least equal to the observed load peak threshold and were often much 
greater. Again, whenever a tentative secondary load peak time slice had 
a predicted load frequency distribution with load values equal to or greater 
than its threshold, it satisfied the peak criteria. Thus, by assuming that 
each predicted load trace was identical in shape to the average load trace 
and by adjusting the load peak threshold for tentative secondary load peaks, 
the predicted load peak criteria were the same as the observed load peak 
criteria. 


b. Calculation of load peak distributions : After the critical time 

slices for a given flight condition were selected, the data for each was 
processed separately. First, each corrected normalized parameter 
distribution was denormalized by multiplying it by the distribution of 
the maximum absolute parameter value. The denormalized distributions 
were found empirically by assuming the independence between each of the 
corrected normalized parameter distributions and the corresponding maxi- 
mum absolute parameter distribution. As a result, the denormalized pa- 
rameter probability distributions expressed the parameter magnitudes 
at each of the critical time slices in terms of the original units. 

For a given maneuver type, the probability of a load peak occurring 
within a specified load range may be calculated by the combination of the 
denormalized peak parameter distributions as indicated by the loads equa- 
tions. The following development illustrates the procedure for the calcu- 
lation of a wing delta load. First, the wing delta load is defined as 

AV W p = c i5( n z - !) + C 17P + C 20P + C 21<1 

Next the probability of a delta load falling in the range from A to B but not 
including B may be expressed as 


Pr [A sAV wp < B ] = Pr [A 5 C 15 (n z - 1) + C 17 p + C 20 p + C 2 x q < B ] 

Then since each denormalized parameter was distributed in 12 ranges and 
each predicted load in 20 ranges, the foregoing equation for AV W p is 
solved for all 20, 736 combinations of the mid-values of the four parame- 
ters in the load equation. Whenever a AV-^-p is outside the threshold value, 
the product of the probabilities of the four parameters An z , p, p, and q is 
calculated. Then the accumulation of these products in the delta load 
ranges gives the probability for each range. Next these probabilities are 
converted to a loads spectrum by multiplying the probabilities by the num- 
ber of maneuvers in the flight condition. As this procedure is repeated 
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for each critical time slice, the number of loads in each range is accumu- 
lated. Finally, since a flight condition defines a steady- state load value, 
the predicted delta loads for all the critical time slices are stored in a 
steady-state range as calculated by = C 15 + C28’ 


C. Data Processing 

As detailed in Appendix G, conventional methods were employed 
to reduce the F-105D eight-channel oscillograph data. This appendix 
also describes the computer programs to form the maneuver model and 
to predict the F-105D load distributions. 


D. Computer Program for Maneuver Pattern Recognition 

How practically the maneuver model technique can be used to pre- 
dict the peak loads on the structural points of an aircraft depends largely 
on the degree to which all phases of the data collecting and processing can 
be automated. Soon to replace oscillograph recorders for large data 
samples, the magnetic tape recorders will markedly advance this auto- 
mation. With such recorders adapted to write data on computer-com- 
patible tapes, the manual aspects of data editing and reading will be 
eliminated. However, before the maneuver model technique can be applied 
in the computer processing of magnetic tape data, the computer must have 
the capability of identifying the maneuver types and of processing their data 
separately. Therefore, to determine the feasibility of automatically iden- 
tifying maneuver types, the current research developed a computer pro- 
gram for maneuver pattern recognition. 

This program was designed to process a tape with data simulating 
that recorded during flight. The data would consist of uncalibrated digi- 
tizations of the eight recorded parameters (n x , ny, n z , p, q, r, airspeed, 
and altitude), each sampled about five times per second. Then the output 
of the program would be a chronological printout of the interesting parame- 
ter values in each maneuver followed by the beginning and ending times and 
the type of the maneuver. The maneuver pattern definitions used in devel- 
oping the computer program and a simplified flow chart are included in 
Appendix D. 
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SECTION II 


RESULTS 


A. Normalized Data by Maneuver Type 

For each of the 23 maneuver types observed in the F-105D data, 
average normalized time histories were prepared for the nine parameters: 
n x , n An z , p, q, r, p, q, and r. As discussed in Section I.B.4.b, these 
time histories were adjusted to yield corrected average normalized pa- 
rameter time histories. For each of the 23 maneuver types, Figure E-l 
in Appendix E pairs the uncorrected and corrected plots for each of the 
nine normalized parameters. 

The number of normalized parameter distributions about each of the 
average normalized values was too large (23 maneuver types x 9 parame- 
ters x 25 time slices = 5175) to permit inclusion in this report. To illus- 
trate such distributions, however, Figure E-2 in Appendix E for the de- 
scending left turn maneuvers shows the corrected normalized distributions 
for the nine parameters at the four critical time slices. As seen here, the 
time slices had average load peaks as follows: the first (No. 5) in both the 
horizontal and the vertical tail load; the second (No. 12) in both the wing 
and the fuselage load; the third (No. 17) in the horizontal tail load, and the 
fourth (No. 19) in the vertical tail load. 

Some statistical tests were made to ascertain whether certain nor- 
malized distributions could be combined. If these tests showed equality 
of most of the normalized parameter distributions between two or more 
maneuver types, the normalized data could be combined into a larger data 
sample, and such results would suggest also testing the maximum absolute 
parameter distributions for equality. If the latter distributions also tested 
equal for the same maneuver types, the number of maneuver types could 
be reduced in the statistical model. 

i 

Since turns provide the largest number of possible maneuver type 
combinations, they were selected for the initial statistical tests. First, 
descending left, left, and ascending left turns were compared with the cor- 
responding right turns at several critical time slices. Although the pa- 
rameters ny and r had a directional difference,, they could be compared by 
inverting the distributions which were stored in ranges symmetrical about 
the normalized zero. The distributions of the normalized parameters n y , 
An z , q, q, and r were tested for equality by the Kolmogorov-Smirnov test 
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at the 0.05 significance level. The tests on the normalized n-y distribu- 
tions showed none equal, whereas the tests on all the q normalized distri- 
butions could not reject equality of the distributions. From a total of 21 
tests on An z distributions, 8 tested not to be equal. Therefore, the nor- 
malized distributions for left turn maneuver types should not be combined 
with those for the corresponding right turn maneuver types. 

Next, statistical tests were made to compare the descending left, left, 
and ascending left turns with each other and, similarly, the descending 
right, right, and ascending right turns with each other to judge whether the 
altitude change could be removed from the criteria to select turn maneuver 
types. Again the tests rejected equality of some, but not all, of the nor- 
malized parameter distributions. Moreover, since no two turn maneuver 
types in the data sample had most of their parameter distributions testing 
equal, none of the maneuver types should be eliminated in future studies. 


B. Maximum Absolute Parameter Distributions 

As obtained from the normalizing of the F-105D maneuver data, the 
maximum absolute parameter distributions are the normalizing factor dis- 
tributions used to denormalize the normalized data preparatory to pre- 
dicting the F-105D load peaks. In addition to the breakdown according to 
maneuver type, the maximum absolute parameter distributions were 
classified by mission segment and flight conditions (combinations of gross 
weight and Mach number) which most affected the coefficients in the loads 
equations. A total of 368 combinations of maneuver type, mission seg- 
ment, and flight condition were used to predict the F-105D loads distri- 
butions, each combination requiring a set of nine maximum absolute 
parameter distributions. 

In the selection of the flight conditions, a preliminary survey of the 
loads equations indicated that gross weight was the most important parame- 
ter and that Mach numbers between 0.9 and 1.2 caused a severe change in 
the horizontal tail loads equation. Because most of the practice external 
stores on the F-105D's were small and light, stores configuration was not 
an important parameter for the F-105D peacetime data. Altitude changes 
also did not seriously affect the loads coefficients. Four mission segment 
groups were formed as follows: (1) ascent, cruise-out and refueling; (2) 

low -angle bombing, high-angle bombing, GAM delivery, special weapons 
delivery, rockets, and air-to-ground gunnery; (3) loiter, air-to-air 
gunnery, and training; and (4) clean cruise and descent. Ranges of gross 
weight and Mach number were grouped into flight conditions so that the 
number of maneuvers for each combination would be as large as possible 
and never less than three. 
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To illustrate the general level of maximum absolute parameter 
values, Table F-l in Appendix F, a computer printout, lists the composite 
distributions for each of the nine parameters versus maneuver type. The 
format limitation of the computer printout required printing each of the 
distributions in two tables: one for maneuver types 1 through 20, and the 

second for the remaining maneuver types. The symbols D/S and D/S/S 
denote degrees per second and degrees per second per second. Section 
I. B.3 identifies the codes for the maneuver types. The number of ma- 
neuvers in the 450 hours of F-105D data totaled 12, 873. 


C. Comparison of Predicted and Observed Loads 

The predicted and observed loads are compared in Figures 4 through 
7. These figures show composite curves of the number of load peaks per 
thousand flight hours above each level of maneuver load for each corre- 
sponding range of level flight load. The predicted and observed curves 
match quite closely for fuselage loads (Figure 4), wing loads (Figure 5), 
and horizontal tail loads (Figure 6). However, the appreciably lower pre- 
dicted curve in Figure 7 for the vertical tail loads indicated that the pre- 
diction technique underestimates the vertical tail load distribution at each 
level . 


To find the source of this underestimation, the predicted and ob- 
served vertical tail load peaks for each maneuver type were examined. 

The findings revealed that the yaw and wing rock maneuvers contributed 
most of the large vertical tail load peaks but that the percentage of error 
in these peaks was about the same as that in the peaks of the other maneu- 
ver types. However, it was also found that the number of vertical tail load 
cycles among the observed yaw and wing rock maneuvers varied from two 
to four per maneuver. On the other hand, the normalized data never pre- 
dicted more than two and a half such cycles per maneuver for these ma- 
neuver types. Apparently, the manual editing should have used a more 
stringent definition of the number of load cycles per maneuver for those 
maneuver types. A slight change in this definition should increase the 
accuracy of the predicted vertical tail load peaks by more than 50 percent 
in the higher load ranges. 

Tables G-l through G-4 in Appendix G are computer tabulations of 
observed and predicted distributions of fuselage, wing, horizontal tail, and 
vertical tail load peaks for each of the 23 maneuver types. However, be- 
cause of the relatively rare occurrence of maneuver types 20, 21, 24, and 
25, their distributions were combined in the tabular presentation. Each 
load heading in these tables represents the lower limit of a range. 
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NUMBER OF PEAKS ^ AVF PER IOOO HOURS 



-2 -6 -10 -14 -18 -22 -26 -30 


FUSELAGE SHEAR MANEUVER LOAD (AVF) IN THOUSANDS OF POUNDS 

Figure 4. —Plots of composite predicted and observed fuselage 
load peaks per 1000 hours 
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NUMBER OF PEAKS* AVW PER IOOO HOURS 



WING SHEAR MANEUVER LOAD (AVW) IN THOUSANDS OF POUNDS 


Figure 5. —Plots of composite predicted and observed wing load 
peaks per 1000 hours 
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Figure 6. —Plots of composite predicted and observed horizontal 
tail load peaks per 1000 hours 




NUMBER OF PEAKS* WT PER IOOO HOURS 



VERTICAL TAIL SHEAR MANEUVER LOAO (WO IN HUNDREDS OF POUNDS 

Figure 7. —Plots of composite predicted and observed vertical tail 
load peaks per 1000 hours 
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D. Pattern Recognition Evaluation 


Three tapes with continuous uncalibrated digital data of .about five 
samples per second for each of the eight recorded parameters (n x , ny, 
n z> P> r, airspeed, and altitude) were used to evaluate the computer 
program for pattern recognition. These tapes contained actual maneuver 
data from fourteen F-105D flights with sections of interpolated data in- 
serted between the maneuvers to simulate the continuous data normally- 
recorded during a flight. As a norm for the comparison of the computer 
recognition, an independent set of pattern recognitions was acquired by 
manually editing the oscillograms corresponding to the tape data. After 
314 maneuvers of various types were recognized by the manual editing 
and 277 by the computer processing, the numbers of individual types were 
grouped according to similar maneuver types; for example, the numbers 
for the six individual types of turn maneuvers were summed under "turns. " 
Tables 7 and 8 summarize the evaluation. Before interpreting the two 
tables, note the apparent discrepancy between the listing of 284 computer 
recognitions (254 correct plus 30 incorrect) in Table 7 and the listing of 
277 in Table 8. This discrepancy was due to several instances where the 
computer recognized two maneuvers as one. With this understanding, 

Table 7 shows the correlation of the computer recognitions with the manual 
recognitions. As seen here, 254 of the computer recognitions were cor- 
rect, 30 were incorrect, and another 30 were not recognized; in other 
words, 81 percent of the 314 manual recognitions had a correct counter- 
part among the computer recognitions, 9. 5 percent had an incorrect cor- 
respondence, and another 9. 5 percent did not sufficiently conform to the 
defined maneuver patterns to permit computer recognition of the maneu- 
ver type. Table 8 summarizes the validity of the 277 computer recog- 
nitions. As stated here, 254 or 92 percent of these recognitions were 
correct. In addition, the computer separated 91 other short sections of 
data but did not recognize any maneuver types in them. For each of these 
sections, however, the computer printed comments identifying those com- 
binations of parameter patterns that made computer recognition impossible 
These 91 sections contained the 30 patterns listed in Table 7 as not rec- 
ognized by the computer. 

In summary, excluding the 30 patterns not recognized by the com- 
puter, or reducing the 314 manual recognitions by this figure to 284, the 
computer correctly recognized 90 percent or 254 of the patterns correctly. 
Now to supplement the computer recognitions by the manual review of the 
printout, or to recognize the remaining 30 patterns (less than 10 percent 
of the total), would require much less than 10 percent of the time expended 
in a completely manual recognition since the computer printed the above- 
mentioned comments for each of the 91 data sections in which the missing 
patterns were contained. 
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TABLE 7 


COMPARISON OF MANEUVER PATTERNS RECOGNIZED BY COMPUTER 
EDITING WITH THOSE RECOGNIZED BY MANUAL EDITING 


Maneuver 

Manual 


Computer Recognitions 

Type 

Recognitions 

Correct 

Incorrect Not Recognized 

Turns 

215 

188 

11 

16 

Pull-ups 

50 

37 

11 

2 

Rolls 

7 

3 

3 

1 

Thrust & Drag 

6 

1 

-- 

5 

Cloverleaf 

10 

7 

3 


Wing Rocks 

13 

10 

1 

2 

Yaws 

1 1 

8 


3 

Pitch- Downs 

2 

-- 

1 

1 

Totals 

314 (100%) 

254 (81%) 30 (9- 5%) 

30 (9. 5%) 


TABLE 8 

VALIDITY OF MANEUVER PATTERNS RECOGNIZED 
BY COMPUTER EDITING 


Maneuver 

_ . T yp _- e _ 

Turns 

Pull-ups 

Rolls 

Thrust & Drag 
Cloverleaf 
Wing Rocks 
Y aws 

Pitch -Downs 
Totals 


Total 

Recognitions 

202 

40 

3 

4 
10 
10 

8 


Correct 

Recognitions 

188 

37 

3 

1 

7 

10 

8 


Incorrect 

Recognitions 

14 

3 

3 

3 


23 (8%) 


277*(100%) 254 (92%) 


In addition to treating the data in which these recognitions were made, 
the computer separated 91 other short sections of data but did not 
recognize any maneuver types in them. As indicated in Table 7, the 
manual editing recognized 30 maneuvers in these sections. 
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SECTION III 


CONCLUSIONS AND RECOMMENDATIONS 


On the basis of the successful prediction of maneuver load peak dis- 
tributions on the fuselage, wing, horizontal tail, and vertical tail of the 
F-105D from 450 hours of eight-channel in-flight data collected during 
peacetime operation, it is concluded that — 

(1) The maneuver model loads prediction technique can be 
adapted to predict such load distributions on a large- 
scale data reduction basis. 

(2) Each of the 23 types of maneuvers observed in the 
F-105D data can be represented by an average nor- 
malized time history and a set of 25 normalized 
distributions for each parameter. 

(3) The normalized data and the set of normalizing factors 
can be recombined to calculate (or predict) maneuver 
fatigue load spectra having accuracies consistent with 
a preliminary fatigue load analysis. 

Prepared independently of the loads prediction development, the 
computer program for pattern recognition can automatically recognize 
and classify maneuvers in the digital time histories of eight-channel data. 
The program recognized 90 percent of all maneuvers in fourteen recorded 
flights of the F-105D, and correctly classified the type of 90 percent of 
the maneuvers recognized. 

In the light of the foregoing results, it is recommended that — 

(1) The hypothesis that the normalized data is independent 
of the aircraft type should be tested by using the F-105D 
normalized data to predict loads on another aircraft type. 

(2) The pattern recognition computer program should be 
further developed to permit inputs of various types and 
formats and to yield outputs in a format compatible with 
the loads prediction computer programs. 

(3) The pattern recognition computer program should be 
tested on some actual magnetic tape flight data and on 
data from another aircraft type. 
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APPENDIX A 


LOADS EQUATIONS AND OBSERVED LOAD PEAK DEFINITIONS 


Development of Loads Equations 

Several simplifying assumptions were made to derive a practical 
set of loads equations. Such assumptions would not likely affect the 
accuracy of the resultant computations seriously since the accuracy was 
already limited by the accuracy of the recorded data, the available wind 
tunnel aerodynamic data, and the available inertia data. These assump- 
tions were as follows: 

(1) The airplane responds as a rigid body. 

(2) The airplane center of gravity is fixed. 

(3) All aerodynamic coefficients are linear functions of a, 3> ix> 
6^, and 

(4) Body axis aerodynamic coefficients are equal to the stability 
axis coefficients given in Reference 11. 

(5) The shape (not the magnitude) of the pressure distribution 
and the center of pressure remain fixed on each aerodynamic 
and control surface. 

(6) The airplane mass distribution is always symmetrical about 
the x, z plane, that is I yz = I = 0. 

(7) The spoiler deflection angle is a linear function of the aileron 
deflection angle for each combination of Mach number and 
altitude. 

(8) The aerodynamic forces on the external stores are negligible. 

(9) The airplane is not exposed to gust; that is, it always flies 
in a uniform free stream. 

These assumptions did not compromise the study objective since 
the intent was to determine the feasibility of applying a statistical cal- 
culation technique to flight loads data, not to calculate loads with strict 
accuracy. The only requirements, therefore, were that the equations 
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APPENDIX A.— Continued 

be approximately correct and that the equations used for the "observed" 
loads also be employed for the statistically predicted loads. 

The first step in deriving a set of loads equations requires setting 
up a balance between the total aerodynamic and inertia forces and be- 
tween the total aerodynamic and inertia moments acting on an airplane. 
Since the study excluded forces in the longitudinal direction, the problem 
was limited to five degrees of freedom with the five unknowns CL, 3> i'p > 
6 a , and 6 R . The five force and moment equations derived for these 
unknowns are as follows: 


n v W 


(1) C Y6a - 5 a+ c Y 6r • 6 R+ c Y p P = ■ 76578 - L ~ z - • 18045 C Y r ^— r 


fa 

V, 


n z W 


(2) Cr • a+ C L . • i T = .76578 

U. u ^ 


(3) C • 6 a + C t • 6 R +C^ • 8= [.00038268 I x p - .0000066791 (I y -I z ) qr 


'If, “ A ' K ’ " 2 

6 A °R P Ve 


.00038268 1 r 000006679 1 I xz pq] - • 18045 ^ - r 

r V c 


- . 18045 C. p 

V e 


C • a + C m . • i T = — [. 001 1637 I q 00002031 1 (I., - I x ) pr -.00002031 1 I xz (r 2 -p 2 )| 


141 -m a “ + ‘T- 2 l " 

v e 


-C m - .059340 C m 


So 


"q v e 


(5) C n6 • 6 A +Cn 6R - 6 r + C n ^ • 8 = — 2 [. 00038268 I z r - .0000066791 (I x - I y ) pq 


- .00038268 I xz p + .0000066791 I xz qr] 

- . 18045 C n r - . 18045 C n p 

v e p v e 


The four shear load equations were formulated as follows: 


Fuselage shear load at Station 277 c 

V F . V. 2 [: 02604 CLo ■ a 4 . 0 1302 * . 00444 C^S-,] 

- 6444 n z + 76. 45 q - 1.333 qr 
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APPENDIX A. —Continued 


Vertical tail shear load at Water Line 38: 


V VT = V 2 r 1 . 302 C Y • 6„ + I. 302 C Yo • g f . 2355 C Y ^—r + . 0523 C, ,, 1 

VT e L y 6 r « Y Pvt r v e p v e j 


v/o 


- 304. 7 n y - L 280 p 4- . 02 234 qr A 3 . 166 r 4 .05 52 pq 


Right horizontal tail shear load at Buttock Line 29 : 

V ht = v e 2 L 380 c m„ +.651 C L . (o' • a+i T ) -.0390C. p 
HT e L °T-TT u lT Cl 1 y 


- . 01684 C m q] - 140. 5 n z - . 002 1 03 (p 2 + q 2 ) + 1 . 680 q 02936 rp 
Q y J 


4 . 3602 p 4 . 00629 qr 

Wing shear load, at right wing Station 136.6 : 

V w .= V 2 f. 1511 (Cl ■ a + C, ) - 1 . 428 Ci, • 6 a - . 550 C . • 6 A 

W 6 L a L »WB *6 A ^'6 s A 

- . 1703 ~~ p]- (460 +W 6 ) n z - . 000031502 W 6 (p 2 4 q 2 ) 

+ (2. 01962 + . 0035069 W 6 ) q 

- (. 03524^ + . 000061207 W^) rp + (3. 49564+ . 0076751 W 6 )p 
+ (. 0610105 + . 000133956 W 6 ) qr 
(where W^ = weight on Station 6, right outboard pylon) 

The solutions of the five unknowns from the foregoing force and 
moment equations and then their substitution in the shear loads equations 
yielded the following loads equations: 


Vyrp - Cj^ny + C 2 P + C 3 r 4- C 4 P -f- C 5 r + C^qr + Cypq 


Vp = Cgn z + C 9 q + Cioq + Ciipr + C 12 qr + C^r 2 - p 2 ) + C 


14 


V W “ ^lb n 7 . + c l6 n y + ^ 17 P + C igq + C19J: + C 2 qp + C 2 iq + C 22 r 
+ C 23 pr + C 24 qr + C 25 (p 2 + q 2 ) + C 2 g(r 2 - p 2 ) + C 27 pq + C 2 g 
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V HT = C 29 n z + C 30 p + c 31 q + c 32 p + C 36 q + C 37 rp + C 38 qr 
+ c 39 (p 2 + q 2 ) + C 40 (r 2 - p 2 ) + C 41 

In these equations the coefficients C} through C 41 are functions of 
the aerodynamic coefficients, Mach number, altitude, dynamic pressure, 
weight, moments of inertia, and so on. The aerodynamic coefficients 
used in the loads equations were based on wind tunnel data given in Ref- 
erence 11. Since this source did not contain sufficient information to 
convert the stability- axis coefficients to body-axis coefficients, the 
stability-axis data was used to approximate the aerodynamic forces and 
moments in the body axis system. As stated above, the inertia and 
weight data was, for the most part, estimated from information in Ref- 
erence 8 and from the physical dimensions of the airplane in Reference 9. 

From the parameter measurements, the computer calculated at 
1 / 5 -second intervals values of Vp, Vyj, Vjj'p* and V^r. These load 
values are hereafter referred to as '’observed loads. " 


Definitions for Observed Load Peaks 

It was noted that the "observed load" deviated about a steady load 
value corresponding to straight and level (n z = 1) flight. By setting 
n z = 1 and n y =p = q= r = p = q= r= Oin the loads equations, the fol- 
lowing steady load equations were obtained: 

V VT v = 0 

VT (n z = 1) 

V F (nz = 1) = C 8 + C 14 

V W(n z= i) = C 15 + C 28 

V HT (nz = i) = C 29 + C 41 


The "observed load" deviations, or delta loads, were then defined as 


AV vt - V VT - v VT( nz = 1) ' V VT 


AV F = V F 


V 


AV W = V W 


F (n z = 1) 
(n z = 1 ) 


v w 


av ht = v 


HT - V HT (nz = i) 
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Given these equations, an "observed load" peak could then be de- 
fined as a delta load beyond a preset threshold and with a rise and fall 
(or fall and rise) equal to or greater than both the threshold value and 
50 percent of the peak delta load. The following lists the thresholds 
which were defined as about 10 percent of the design loads in Reference 12. 

-6000 < AVp < + 6000 lb. 

-1800 < AV vt < + 1800 lb. 

-800 < A ht < 800 lb. 

-4000 < AV W < + 4000 lb. 

On the basis of these thresholds, the computer tabulated the delta 
and steady loads of each "observed load" peak for later comparison with 
the statistically predicted load peaks. 

Because the statistical prediction of loads involves predicting the 
probabilities for all possible parameter combinations, the number of 
independent parameters in each load equation should be reduced to that 
minimum which would yield an accuracy consistent with that expected in 
the load prediction technique. Consequently, on the basis of the relative 
magnitude of these terms, the equations to predict loads were reduced 
to the following approximations, where the subscript "P" denotes an 
approximated expression: 

V VTP = C l n y + C 2 P + ^5 r + C 7 Pq 
V FP = c 8 n z + Ci 0 q + 

V WP = c 15 n z + C 17 P + C 2 qP + C 2 iq + C 28 
V HTP = C 29 n z + C 3 iq + c 36q + c 37 r P + C 41 
AV VTP = V VTP 
AV FP = V FP ' V F(n z = 1) 

AV WP = V WP ‘ V W(n z = 1 ) 

AV HTP = V HTP “ V HT(n z = 1 ) 
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APPENDIX B 

PARAMETER PATTERNS FOR BASIC MANEUVERS 


The application of the maneuver model to predict structural load dis- 
tributions requires that maneuver types be recognizable in eight-channel 
data. The parameters most indicative of maneuver types are the three 
angular velocities, p, q, r; normal acceleration, n z ; lateraL acceleration, 1 
n yj and altitude. The airspeed and the longitudinal acceleration, n x , can 
sometimes support these parameters. The trace patterns indicative of the 
basic maneuver types are described below. 

The turn maneuver may be identified by the combination of the follow- 
ing trace patterns: (1) a long positive peak in the normal acceleration, n z , 

trace; (2) depending on the turn going either right or left, a long positive or 
negative peak in the yaw rate, r, trace; (3) again depending on the direction 
of the turn, an early positive or negative peak followed by a late peak of 
opposite sign in the roll rate, p, trace. In addition, a long positive peak 
similar to that in the normal acceleration trace appears in the pitch rate, 
q, trace. The trend in the altitude trace indicates that the turn is ascend- 
ing, descending, or level. Figures B-l and B-2 show oscillograph record- 
ings of a descending left turn and a descending right turn, respectively. 

As variants of the turn maneuver, the left and right cloverleaf maneuvers 
have the following trace characteristics: the altitude increases and then 

decreases. However, rather than one positive normal acceleration peak 
as in the turn maneuver, these maneuvers have two such peaks and the high 
point in the altitude trace corresponds with the dip between them. The roll, 
yaw, and pitch rate traces have the same patterns as in a turn maneuver. 

A pull-up maneuver may be classified as either a rolling or a sym- 
metrical pull-up depending on whether or not it includes a roll rate peak. 

A symmetrical pull-up maneuver may be identified by the combination of 
the following trace patterns: (1) a large positive peak in the normal accel- 

eration trace; (2) a large positive peak in the pitch rate trace occurring 
simultaneously with the former peak; and (3) an increasing rate of climb or 
slope in the altitude trace. Figure B-3 shows an oscillograph recording of a 
symmetrical pull-up. Besides all the characteristic trace patterns of the 
symmetrical pull-up, a rolling pull-up maneuver has a large positive or 
negative roll rate and yaw rate peak, both occurring during the duration of 
a large normal acceleration peak. The sign of the roll rate and yaw rate 
peaks indicates the direction of the roll. Figure B-4 shows an oscillograph 
recording of a right rolling pull-up. Although the rolling pull-up maneuver 
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is quite similar to a symmetrical pull-up maneuver followed by an ascend- 
ing turn maneuver, the two pull-up maneuver types can be distinguished by 
noting when the roll rate peak begins. If it begins while the normal accel- 
eration is still high, the maneuver is a rolling pull-up. But, if the peak 
begins when the normal acceleration has returned close to a 1.0 value, 
there are two maneuvers, that is, a symmetrical pull-up and a turn. As 
a variant of the symmetrical pull-up maneuver, the symmetrical pitch- 
down maneuver has a decreasing altitude and a negative normal acceler- 
ation. As another variant of the symmetrical pull-up maneuver, the 
inside-loop maneuver has a sustained pitch rate deflection until the air- 
plane completes the loop. During the inside-loop maneuver, the altitude 
increases and then returns close to the altitude at the start of the maneu- 
ver. Also during this maneuver, the normal acceleration has an early 
positive peak, a small negative peak at the inverted position,, and a late 
positive peak. 



Bottom Reference 


Figure B-l. —Oscillogram showing descending left turn 
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START 
PULL - UP 


START 

TURN 


Figure B-4. —Oscillogram showing right rolling pull-up 


The yaw maneuver is characterized by a deflection in the yaw rate 
trace and a large deflection in the lateral acceleration trace. None of the 
other parameters vary significantly. Figure B-5 illustrates a yaw maneu- 
ver. During the recording of the F-105D data, a yaw maneuver was often 
performed early in a flight by producing a right and left yaw in quick suc- 
cession (that is, a "rudder kick") to test the rudder control system. 

The acceleration and deceleration maneuvers indicate an abrupt 
power change or the use of either an afterburner or a dive-brake system. 

A rapid increase or decrease of the longitudinal acceleration characterizes 
these maneuvers. These changes are of relatively short duration and end 
as the longitudinal acceleration returns to a normal value. While the air- 
speed trace increases or decreases, none of the other parameters vary 
appreciably. Figure B-6 shows a deceleration maneuver. 
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APPENDIX B. —Continued 

The roll maneuver is characterized by a long peak in the roll rate 
trace, p. The yaw rate trace, r, moves first in one direction and then in 
the other because of the induced yaw and the pilot's subsequent action to 
correct for the induced yaw. The normal acceleration trace has no sig- 
nificant activity. Figure B-7 shows a roll maneuver. As seen here, a 
180-degree roll followed a 180-degree pitching maneuver which ended 
with the aircraft inverted. Related to the roll maneuver is the wing rock 
maneuver which is normally characterized by several cycles of roll rate 
oscillations with roll angles less than 90 degrees. In the wing rock ma- 
neuver, the yaw rate also has oscillations which slightly lag those of the 
roll rate. Also in the wing rock maneuver, the lateral acceleration re- 
sponds inversely to the slope of the yaw rate trace; that is, the lateral 
acceleration is positive when the yaw acceleration is negative and the 
lateral acceleration is negative when the yaw acceleration is positive. As 
variants of the basic roll maneuvers, the right and left four-point roll 
maneuvers each have four distinct peaks in the roll rate as it remains posi- 
tive or negative throughout the maneuver. Like the wing rock maneuvers, 
the four-point roll maneuvers have the roll and yaw rates oscillating in 
several cycles with the yaw rate oscillations slightly lagging those of the 
roll rate and the lateral acceleration responding similarly to the yaw rate 
deflections . 





Figure B-7. —Oscillogram showing roll 
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APPENDIX B. —Concluded 

Although the maneuver types described above include all those ob- 
served in the available F-105D data, maneuvers not covered by these types 
may be expected in other data. Then, either the description of an existing 
type will have to be enlarged or a new maneuver type defined. 
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APPENDIX C 


DATA PROCESSING 


Data Reduction Methods 

The following paragraphs describe the major steps to reduce the 
F-105D flight loads data used in this study: editing, reading, quality 

control checking, and computer output checking. 

1 . Editing . — At the outset of the editing, each oscillogram was re- 
viewed to detect any evidence of instrument malfunctioning and to check 
the adequacy of the calibration signals. In addition, the supplemental data 
sheets were inspected to confirm the completeness and correctness of 
their information. 

Next the oscillograms were scanned to mark off the sections to be 
read. Except for the few instances where turbulence obviously caused the 
trace deflections, such sections consisted of all maneuvers where any pa- 
rameter exceeded its normalizing threshold (as defined in Section I.B.4.a) 
and those maneuvers where all parameters were within the thresholds but 
where the pattern of a maneuver type was recognized. Then all selected 
sections were identified by a base, mission type, mission segment, ma- 
neuver type, takeoff and landing configuration, takeoff and landing gross 
weight, takeoff and landing fuel weight, store weights, takeoff barometric 
pressure, and calibration information. 

2. Reading . — Two procedures governed the reading of the selected 

maneuver sections: one for sections where at least one parameter ex- 

ceeded its normalizing threshold, and the second for sections where no 
parameter exceeded its normalizing threshold. According to the first 
procedure, the n x , n y , n z , p, q, and r traces were digitized at 0.2- 
second intervals, and the static and dynamic pressure traces at 2.0- 
second intervals. The computer later processed this data to calculate 
the time histories of the observed loads and then to form the time his- 
tories of the normalized parameters. 

According to the second procedure, only three points of each pa- 
rameter were digitized in each maneuver section: the first, the maximum, 
and the last. The computer later processed this data to determine the 
total number of maneuvers of each type and to complete the maximum ab- 
solute parameter distributions at the lower levels. Although the computer 
calculated the observed loads for these maneuvers, none were large enough 
to be classified as observed load peaks. 
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Finally, a preliminary printout of the digitized data was reviewed 
to detect any obvious reading errors. 

3. Quality control checking. — The quality control check served to 
hold within acceptable limits the inevitable small reading errors resulting 
from thousands of measurements taken to the nearest 0.005 of an inch. 

For this check, therefore, randomly selected digitized readings from each 
flight were compared with the manual measurements of the corresponding 
points on the oscillogram. Then the differences (or errors) were used to 
establish and then maintain acceptable limits for the reading deviations. 
Whenever the data of an entire flight was found acceptable, it was for- 
warded for computer processing; otherwise, it was reread. 

For each of the eight digitized parameters. Table C-l lists the mean 
and the standard deviation (o ) of the reading errors found in the quality 
control check of data from 378 flights. On the assumption of a Gaussian 
distribution of the reading errors, 95 percent of all measurements should 
not differ by more than 2a from the correct value. 


TABLE C-l 

QUALITY CONTROL READING ERROR STATISTICS 

Standard Deviation O 


Parameter 

Mean of Reading Error 

of Reading Error 

n x 

0. 0003 g 

0. 006 g 

n y 

-0. 0006 g 

0. 006 g 

n z 

-0. 004 g 

0. 04 g 

P 

-0. 08°/ sec 

0. 9°/ sec 

q 

-0. 014°/sec 

0. 14°/sec 

r 

-0. 014°/sec 

0. 1 5°/ sec 

Altitude 

4 ft. * 

84 ft. * 

Airspeed 

0. 1 knot** 

2. 3 knots* 


* at 5000 feet 

** at 400 knots 
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4. Computer output checking . — The computer program was de- 
signed to print out comments indicating any data having unusual parameter 
values or trace patterns when compared with the F-105D theoretical flight 
envelope and with those expected for a particular combination of maneuver 
type, mission type, mission segment, and configuration. In the manual 
review of the printout, the parameter values and supplemental data were 
checked to detect all large reading errors and to verify the proper classi- 
fication of all maneuvers. The data of each flight was reprocessed until 
the printout showed no errors. 


Computer Programs to Form Maneuver Model and 
to Predict F-105D Load Distributions 

According to the definitions and procedures described above, the 
first computer program (Phase I) processed the digitized parameter data 
and supplemental data to yield the following separate distributions: (1) 

normalized parameter values (as defined in Section I. B. 4. a) by maneuver 
type, (2) maximum absolute parameter values by maneuver type and flight 
condition, and (3) observed load peak values (as defined in Appendix A) 
by maneuver type and flight condition. The results of Phase I were printed 
for the computer output check and transcribed on magnetic tape for further 
processing. 

The second computer program (Phase II) combined all acceptable 
flight data on the Phase I output tapes onto a single Phase II output tape 
to serve as the input for the third and fourth computer programs. 

The third computer program (Phase III) printed tabular distributions 
of the data contained in the Phase II output tape. 

The fourth program (Phase IV), the loads prediction program, first 
sorted and stored the proper normalized distributions, the maximum abso- 
lute parameter distributions, and the observed load distributions. A card 
input established the maneuver types to be predicted and the exact ranges 
for the flight condition variables defining the flight condition. Next the 
fourth program calculated the constants for the loads equations and the 
average load time histories, as described in Section I. B. 5. a. Then it cal- 
culated the predicted load probabilities and multiplied them by the number 
of maneuvers in the given flight condition to accumulate load frequencies. 

The fourth program had three levels of output. Level 1 gave the 
predictions for one flight condition and the associated values calculated 
for each of the flight condition variables. Levels 2 and 3 provided com- 
posite tables of predicted loads. Level 2 summarized the loads from a 
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specific mission segment group, and Level 3 summarized all Level 2 
tables giving predicted loads for all flight conditions of a maneuver type. 

At each of the three levels, the tables for observed loads were printed for 
ready comparison. Also available on option at Levels 2 and 3 were com- 
posite probability plots presenting the cumulative probability of loads for 
both the predicted and the observed peaks. 

Besides the above-mentioned tape input for the fourth program, a 
card-input capability was provided to transfer statistical distributions 
directly from cards and to store them in the computer memory. This 
added capability gave a greater flexibility in studying prediction techniques. 
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PATTERN RECOGNITION DEFINITIONS 


The major effort in developing a pattern recognition computer pro- 
gram is to rigorously define the characteristic parameter patterns in the 
time history of each maneuver type. To effect the automating of the 
manual editing as much as possible, all maneuvers should be classified. 

In addition, the parameter pattern definitions should preclude misclassi- 
fications since such would reduce the validity of the entire data sample 
and, therefore, be worse than no classification at all. Consequently, 
the objective of the maneuver pattern definitions was twofold; (1) the 
classification of as many maneuvers as possible in the data, and (2) the 
rigorous definition of parameter patterns for each maneuver type which 
would virtually obviate misclassifications . 

The loads prediction phase of the current study revealed that no 
observed nor predicted load peaks could be found in maneuvers where 
all parameters were below some threshold. In addition, the relatively 
significant but random deflections caused by turbulence often obscured 
and confused the patterns of such maneuvers. Therefore, the maneuver 
pattern definition included the following parameter thresholds which en- 
sured the inclusion of all maneuvers yielding significant load peaks but 
the exclusion of those not yielding such peaks: 

-0. 2 g < n x < 0. 2 g -30°/sec < p < 30°/sec 

-0. 1 g < ny < 0. 1 g -5°/sec < q < 5°/sec 

-1.0 g < An z < 1. 0 g -5°/sec < r < 5°/sec 

In the following presentation of the maneuver pattern definitions, 
the maneuver types are listed generally in the same order that the com- 
puter program evaluated the data sections for maneuver types. As shown 
in the diagram of Figure D-l, the logic flow of the pattern recognition 
program indicates the actual order in which the maneuver types were 
tested in the evaluation of each data section. 

Wing rock maneuver . — This maneuver type must have three roll 
rate peaks in alternately opposite directions. A roll rate peak is defined 
as three consecutive roll rate readings each above 24 counts (approx. 

10° /sec.). The last roll rate peak must occur within 12 seconds (60 
readings) after the start of the maneuver. There must be no sustained 
An z deflection defined as five consecutive An z readings each above 49 
counts (approx. 1.0 g). This maneuver type ends 0.4 seconds after the 
third roll rate deflection falls below threshold or below half of the third 
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Figure D-l. — Logic flow of pattern recognition program 
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peak value, whichever of the two events occurs later. The maneuver is 
classified as a left wing rock if the first roll rate is negative or as a 
right wing rock if the first roll rate is positive. 

Clover leaf maneuver. — This maneuver type must have two sus- 
tained An z deflections. Each is defined as 65 consecutive readings each 
above 14 counts (approx. 0. 3 g). But neither can be an extremely large 
An z defined as five consecutive readings each above 99 counts (approx. 

2. 3 g) since such would more likely reflect a pull-up maneuver. The 
clover leaf maneuver must also have a sustained yaw rate deflection 
with 100 consecutive readings each above 14 counts (approx. 1. 5°/sec. ). 

If the sustained yaw rate deflects negatively, the maneuver is classified 
a left clover leaf; but if it deflects positively, the maneuver is classi- 
fied a right clover leaf. This maneuver type ends 0.4 seconds after 
whichever of the following events occurs last: the second An z deflection 

falls below threshold or below half of the second An z peak value, the 
pitch rate deflection falls below threshold or below half the pitch rate 
peak value, or the yaw rate deflection falls below threshold or below 
half the yaw rate peak value. 

If none of the foregoing patterns have been detected in a data sec- 
tion, the computer program reviews the An z deflections. If five con- 
secutive An z 's occur outside a threshold of i 19 counts (approx. — 0. 45 g), 
the program checks for a pitch-down maneuver when the n z deflection 
is negative or for a four-point roll, a pull-up, or a turn maneuver when 
the deflection is positive. If such An z 's do not appear, the computer 
checks for a roll, a yaw, a deceleration, or an acceleration maneuver. 

Pitch-down maneuver . — This maneuver type must have five con- 
secutive An z readings each below -19 counts (approx. -0.45 g) and five 
consecutive pitch rate readings each below -14 counts (approx. 1. 0°/sec. ). 
In this maneuver type, the pressure altitude normally decreases by at 
least -15 counts (an approx. 700-foot decrease at 5000 feet); if not, the 
program still accepts the pattern as a pitch-down maneuver but adds a 
comment to the classification. This maneuver ends 0.4 seconds after 
whichever of the following events occurs last: the An z deflection rises 

above threshold or above half the negative An z peak value or the pitch 
rate deflection rises above threshold or above half the negative pitch 
rate peak value. 

Roll maneuver . — This maneuver type must have five consecutive 
roll rate readings each above i 74 counts (approx. — 27°/sec. ) and the 
integral of the roll rate trace must be greater than i 100°. Negative 
roll rates indicate a left roll, and positive roll rates a right roll. This 
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maneuver type ends 0. 4 seconds after whichever of the following events 
occurs last: the roll rate deflection falls below threshold or below half 

the roll rate peak value or the An z falls below threshold. 

Yaw maneuver. — This maneuver type must have two opposite yaw 
rate peaks where the first is defined as two consecutive readings each 
abovei 24 counts (approx. — 2. 5°/ sec. ) and the second peak as one reading 
above i 24 counts. There must not be more than 5 seconds between the 
peaks nor more than 12 seconds between the maneuver start and the 
second peak. If the first yaw rate peak is negative, the maneuver is 
a left yaw; if positive, it is a right yaw. This maneuver type ends 0. 4 
seconds after the second deflection returns to threshold or to half the 
second peak value, whichever is later. 

Deceleration maneuver . — This maneuver type must have five con- 
secutive negative n x readings each below -15 counts (approx. -. 05 g) 
and each preceded and followed by at least one positive n x reading more 
than 32 counts (approx. + . 1 g) above the negative n x peak value. In 
this maneuver type, the airspeed normally decreases by -15 counts 
(an approx. 20-knot decrease at 350 knots) or the altitude increases 20 
counts (an approx. 1000-foot increase at 5000 feet); if not, the program 
still accepts the pattern as a deceleration maneuver but adds a comment 
to the classification. This maneuver type ends when the n x deflection 
rises above threshold. 

Acceleration maneuver. — Except for the reversal of the deflection 
signs, this maneuver type has the same pattern as that defined for the 
deceleration maneuver. 

Four-point roll maneuver . — This maneuver type must have four 
roll rate peaks each in the same direction and each with five consecutive 
roll rate readings abovei 98 counts (approx, i 3 5°/ sec. ). Between each 
two adjacent roll rate peaks, the roll rate must fall below 99 counts. In 
addition, this maneuver type must have five consecutive An z readings 
each above 19 counts (approx. 0. 45 g) and five consecutive pitch rate 
readings each above 14 counts (approx. 1. 0°/sec. ). This maneuver 
type ends 0. 4 seconds after the fourth roll rate deflection or the An z 
value falls below threshold, whichever occurs last. Negative roll rates 
indicate a left four-point roll, and positive roll rates a right four-point 
roll. 


At this stage, the computer program determines whether the roll 
rate deflection is sufficient to classify the maneuver a turn or a rolling 
pull-up. If the maneuver has five or more consecutive An z readings 
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each above 99 counts (approx. 2. 3 g), it must also have five roll rate 
readings each above i 59 counts (approx. — 22°/sec. ) to be classified 
either of the two types; but if these An z 's are lower, the roll rates need 
only exceed i 14 counts (approx, i 5. 5°/sec. ). 

Symmetrical pull-up maneuver. —This maneuver type must have 
five consecutive An z readings each above 19 counts (approx. 0.45 g) 
and five consecutive pitch rate readings each above 14 counts (approx. 
1.0°/sec. ). In this maneuver type, no five consecutive roll rate read- 
ings can be outside threshold (i 59 counts if the five or more An z 's are 
above 99 counts or i. 14 counts if the An z ’s are lower), nor can five 
consecutive yaw rate readings be outside threshold (ft 46 counts if the 
five or more An z ’s are above 99 counts or ft 14 counts if the An z ’s are 
lower). The peak An z reading must be at least 25 counts above the 
first An z reading. Finally, the altitude normally has a 15~.count in- 
crease after the beginning of the An z peak; if not, the computer pro- 
gram adds a comment to the symmetrical pullup classification. This 
maneuver type ends 0. 4 seconds after the An z trace falls below thresh- 
old or below half the An z peak or the pitch rate trace falls below thresh- 
old or below the pitch rate peak value, whichever occurs last. 

Right pull-up maneuver . — This maneuver type must have five con- 
secutive An z readings each above 19 counts (approx. 0.45 g) and five 
consecutive pitch rate readings each above 14 counts (approx. 1.0°/sec. ). 
In addition, it must have five consecutive roll rate readings each above 
threshold (i 59 counts if five or more consecutive An z readings are each 
above 99 counts or — 14 counts if the An z 's are lower). The first roll 
rate reading outside threshold must be delayed at least 0. 2 seconds after 
the first An z reading outside threshold. Moreover, this maneuver type 
must also have five consecutive yaw rate readings above threshold (ft 45 
counts if the roll rate is below threshold and the five or more consecu- 
tive An z 's are each above 99 counts; or ft 14 counts if the An z , s are lower). 
Like the first roll rate reading, the first yaw rate reading outside thresh- 
old must be delayed at least 0. 2 seconds after the first An z reading out- 
side threshold. The peak An z value must be at least 25 counts above the 
first An z reading. Finally, the altitude trace normally had a 15-count 
increase after the beginning of the An z peak. If a maneuver without this 
altitude increase has five consecutive An z readings each above 99 counts 
(approx. 2. 3 g), the computer program still classifies it a right pull-up 
and adds a comment. But if a maneuver has neither the altitude increase 
nor these An z , s, the program classifies it a turn and adds a comment. 

The right pull-up maneuver ends 0. 4 seconds after whichever of the fol- 
lowing events occurs last: the An z , the pitch rate, the yaw rate, or the 

negative peak of the roll rate falls below threshold or to half their re- 
spective maximum deflection values. The negative roll rate peak used 
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in this definition of the maneuver termination must be within 25 readings 
before or after the time when the An z and yaw return to threshold. 

Left pull-up maneuver . — Except for the reversal of the roll and 
yaw rate deflections, this maneuver type has the same pattern as that 
defined for the right pull-up maneuver. 

Right turn maneuver. — This maneuver type must have five con- 
secutive An z readings each above 19 counts (approx. 0.45g), five con- 
secutive pitch rate readings each above 14 counts (approx. 1. 0°/sec. ), 
and five consecutive positive roll rate readings each above threshold 
(+ 59 counts if five or more consecutive An z readings are above 99 
counts or + 14 if the An z f s are lower). The first roll rate reading 
outside threshold must occur before or at the first An z reading outside 
threshold. In addition, this maneuver type must have five consecutive 
yaw rate readings each above +14 counts (approx. +1. 5°/ sec. ). The 
positive roll rate deflection must be followed by a negative roll rate 
deflection defined by five consecutive readings each below -9 counts 
(approx. -4. 0°/ sec. ). This maneuver type ends 0. 4 seconds after 
whichever of the following events occurs last: the An z , the pitch rate, 

the yaw rate, or the negative peak of the roll rate trace falls below 
threshold or to half their respective maximum deflection values. For 
the negative roll rate peak used in this definition of the maneuver ter- 
mination, the computer program first searches for one within 25 read- 
ings before or after the time when the An z and yaw rate return to thres- 
old. But not finding such, the program uses that negative roll rate peak 
which occurs first in the maneuver pattern. If the difference is -15 
counts or below, the computer program classifies the maneuver a de- 
scending right turn; if between -15 counts and +15 counts, a right turn; 
and if +15 counts or above, an ascending right turn. 

Left turn maneuver . — Except for the reversal of the roll and yaw 
deflections, this maneuver type has the same pattern as that defined for 
the right turn maneuver. 

Inside loop maneuver . — The few maneuvers of this type were 
noted too late in the effort to be included in the pattern recognition com- 
puter program. 

In its present form, the computer program prints in chronological 
order the data for each section examined for a maneuver. As defined 
above, such data sections include only those where at least one param- 
eter has deflections outside threshold. At the end of each printout of a 
data section, it also prints the maneuver type with or without comments 
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and the beginning and ending times of the maneuver or simply comments 
when no maneuver is recognized. Whenever a recognized maneuver has 
an unusual pattern, the program prints a comment to that effect along 
with the maneuver type. With further development, the program will 
likely print the entire data section only when a comment is generated 
because a data section has either no recognizable maneuver type or an 
unusual pattern in a recognized maneuver type. As illustrated in Fig- 
ure D- 1 , the numbers above 9100 enclosed in circles indicate the places 
where the logic not recognizing a maneuver pattern generates a comment 
giving the reason for no pattern recognition. The numbers between 100 
and 1000 also enclosed in circles represent the places where the logic 
finds a data section with the potential for another maneuver type and 
accordingly transfers it to another phase. The numbers 1 to 2 5 again 
encircled indicate the places where the logic recognizes one of the 
twenty-three maneuver types and generates the maneuver classification. 
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NORMALIZED DATA 


Contents: 

Figure E-l. — Average Normalized Time Histories 
of Parameters in All Maneuver Types 

Figure E-2, — Corrected Normalized Parameter 
Distributions in the Descending 
Left Turn Maneuver 
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TIME SLICE NO. 5 



TIME SLICE NO. 12 




TIME SLICE NO 19 



NORMALIZED N x 

(a) n x — longitudinal load factor 

-2. —Corrected normalized parameter distributions 
in the descending left turn maneuver 
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TIME SLICE NO. 5 



NORMALIZED q 


TIME SLICE NO. 12 



NORMALIZED q 


TIME SLICE NO. IT 



NORMALIZED q 



NORMALIZED q 

(e) q— -pitch angular velocity 

2. —Corrected normalized parameter distributions 
in the descending left turn maneuver (continued) 



























APPENDIX F 


MAXIMUM ABSOLUTE PARAMETER DISTRIBUTIONS 


Contents: 

Table F-l. —Maximum Absolute Parameter Distributions 
by Maneuver Type 
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TABLE F-l. —MAXIMUM ABSOLUTE PARAMETER DISTRIBUTIONS BY MANEUVER TYPE 

( c ) An 7 — incremental normal load factor (continued) 
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TABLE F- 1 . —MAXIMUM ABSOLUTE PARAMETER DISTRIBUTIONS BY MANEUVER TYPE 

(e) q — pitch angular velocity (continued) 
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TABLE G-l. —OBSERVED AND PREDICTED FREQUENCIES OF THE FUSELAGE LOADS 

(continued) 

Maneuver Tvne 4 — Descending right turn 


APPENDIX G. —Continued 
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UENCIES OF THE FUSELAGE LOADS 


APPENDIX G. —Continued 
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TABLE G- 1 . —OBSERVED AND PREDICTED FREQUENCIES OF THE FUSELAGE LOADS 

(continued) 

Manuever Type 6 — Ascending right turn 


APPENDIX G. —Continued 
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TABLE G- 1 . —OBSERVED AND PREDICTED FREQUENCIES OF THE FUSELAGE LOADS 

(continued) 

Maneuver Type 7 — Symmetrical pull-up 


APPENDIX G. —Continued 
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APPENDIX G. —Continued 
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TABLE G- 1 . —OBSERVED AND PREDICTED FREQUENCIES OF THE FUSELAGE LOADS 

(continued) 

Maneuver Type 9— Left rolling pull-up 


APPENDIX G. —Continued 
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TABLE G-l. —OBSERVED AND PREDICTED FREQUENCIES OF THE FUSELAGE LOADS 

(continued) 

Maneuver Type 10 — Right roll 


APPENDIX G. —Continued 
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TABLE G-l. OBSERVED AND PREDICTED FREQUENCIES OF THE FUSELAGE LOADS 

(continued) 


APPENDIX G. —Continued 
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APPENDIX G. —Continued 
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TABLE G-3. —OBSERVED AND PREDICTED FREQUENCIES OF THE HORIZONTAL TAIL LOADS 

Maneuver Type 1 3 — Longitudinal acceleration (continued) 
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TABLE G-4. —OBSERVED AND PREDICTED FREQUENCIES OF THE VERTICAL TAIL LOADS 
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Maneuver Type 3— Ascending left turn 
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